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Abstract 

Diffractive photoproduction of n c is an important process to study the effect of Odderon, 
whose existence is still not confirmed in experiment. A detailed interpretation of Odderon in 
QCD, i.e., in terms of gluons is also unclear. Taking charm quarks as heavy quarks, we can 
use NRQCD and take r\ c as a cc bound state. Hence, in the production of i) c a free cc pair is 
first produced and this pair is transformed into n c subsequently. In the forward region of the 
kinematics, the cc pair interacts with initial hadron through exchanges of soft gluons. This 
interaction can be studied with HQET, which provides a systematic expansion in the inverse 
of the c-quark mass m c . We find that the calculation of the S'-matrix element in the forward 
region can be formulated as the problem of solving a wave function of a c-quark propagating 
in a background field of soft gluons. At leading order we find that the differential cross-section 
can be expressed with four functions, which are defined with a twist-3 operator of gluons. 
The effect of exchanging a Odderon can be identified with this operator in our case. We 
discuss our results in detail and compare them with those obtained in previous studies. Our 
results and those from other studies show that the differential cross-section is very small in 
the forward region. We also show that the production through photon exchange is dominant 
in the extremely forward region, hence the effect of Odderon exchange can not be identified 
in this region. For completeness we also give results for diffractive photoproduction of J/^>. 

PACS numbers: 12.38.-t, 12.39.Hg, 13.60.-r, 13.60.Le 

Key Words: Soft gluon, HQET, diffractive photoproduction of n c , HERA experiment. 
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1. Introduction 



Diffractive photoproduction of a pseudoscalar meson is an interesting process because the 
production is related to the postulated object: OdderonpQ, the partner of Pomeron. The Pomeron 
is even under charge conjugation C, while the Odderon has C = — 1. The exchange of Pomeron 
and Odderon is believed to delivery dominant contributions to hadronic cross section at high 
energy. An interesting review about the history of Odderon and relevant references can be found 
in P]. Giving the importance of its existence in theory, experimentally it is still unsuccessful to 
hunt the effect induced by Odderon. Only one indication for existence of Odderon is found in the 
t-dependence of pp- and pp elastic cross-section [3J. Theoretically the effect of Odderon exchange 
in pp- and pp scattering at high energy has been intensively studied, a recent work and useful 
references can be found in jl]. 

It has been suggested that the effect of Odderon may be detected in diffractive pseudoscalar 
meson production for collision of a hadron with a real- or virtual photon^ Id^ I7| I§1 ITU1 lllj. where 
effects of Pomeron are absent. It is difficult to made predictions for the proposed processes by 
starting from QCD directly. Some models are introduced to make predictions for the ep collider 
at HERA. Unfortunately, the experimental study in photoproduction of 7r°|I2J gives a negative 
result. In jS] it is suggested to look for Odderon in productions of two pions. In j^J E3 E] 
diffractive photoproduction of r\ c are studied, predictions may be made more reliably than those 
for production of light meson, because of the following reasons: The structure of rj c is simpler than 
light mesons, the large quark mass m c enables us to use perturbative QCD at certain level and 
the c-quark content in the initial hadron can be neglected. Therefore the production of r] c can 
be imagined as that the initial photon is split into a cc pair, this pair exchanges gluons with the 
initial hadron h and forms the produced r\ c after exchanges. This is illustrated in Fig.l. Effect 
of exchanging gluons can be thought as the effect of exchanging a Odderon. In [OJ E3] one uses 
perturbative QCD to handle the emission of gluons by the cc pair, at leading order only three 
gluons are emitted. The interaction of the three gluons with the initial hadron h is described by a 
impact factor. Using models for the impact factor one can obtain some numerical predictions. In 
[TT] one uses a set of new Odderon states to account the effect of Odderon exchange and numerical 
predictions can also be made. All of these works deliveries a total cross-section for HERA roughly 
at order of 10 2 pb, but the t-dependence is predicted differently in different works. 
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Figure 1: Typical diagram for exchange of many soft gluons for r) c photoproduction 
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In this work we study the problem from another point of view. We make an attempt to answer 
the question if one can express the differential cross-section in the nearly forward direction with 
quantities which are well defined in the framework of QCD. The answer is closely related to a 
QCD interpretation of Odderon in our case. If r] c is produced diffractively and the beam energies 
are large enough, the exchanged gluons are soft. It is questionable to use perturbative QCD for 
these soft gluons. But, the charm quark can be taken as a heavy quark, the emission of these soft 
gluons can be studied with the Heavy Quark Effective Theory(HQET)[.Q], in which a systematical 
expansion in m^ 1 can be made. Also by taking charm quark as a heavy quark one can describe 
r] c with nonrelativistic QCD(NRQCD) [Tl] . in which a systematic expansion in the small velocity 
v c , which is the velocity of a c- or c quark inside a charmonium in its rest frame, can be made. At 
leading order of v c the cc pair after exchanging soft gluons can be taken as on-shell. Using this 
fact, the study of the emission of soft gluons by the c- or c quark can be formulated to solve a wave 
function of the c- or c quark propagating in a background field of soft gluons. The problem of 
light quarks propagating in a background field of soft gluons have been studied in ^3] i n relation 
to hadron-hadron scattering at high energy. In our case these wave functions can be solved by 
expanding them in m~ l . It is interesting to note that at the limit m c — > oo the exchanged Odderon 
here consists of three soft gluons but in a special gauge. 

The approach outlined above has been used to study similar cases with exchanges of soft 
gluons like J/ty — > 7* + soft pions[ini HZj and diffractive photoproduction of J/^fTS]. In [TH] it is 
assumed that results for the S-matrix element can be obtained by taking two-gluon exchange in a 
special gauge. In this work we derive these results for completeness and without the assumption. 
It shows that for the case of J/\& the exchanging Pomeron consists of two gluons in this special 
gauge. 

In this work we are unable to give numerical results in detail, because the differential cross- 
section is expressed with four unknown functions, which are defined with a twist-3 gluonic operator 
in QCD. However, order of magnitude can be estimated in comparison with the case of J/^f. It 
turns out that the differential cross-section is very small as that predicted by model calculations [3 
[HIJ E] • With our result it can be show that the differential cross-section can be non zero in the 
limit t — > 0, in contrast to the predicted in [HIEIIE], where t is the squared momentum transfer 
between light hadrons. The reason for this is that there are nonperturbative effects which can 
cause the helicity flip of the initial hadron. This will be discussed in detail. Since the differential 
cross-section with exchanges of soft gluons is small, the exchange of a photon, instead of exchanging 
many gluons in Fig.l., can have a sizeable contribution. At first look, the amplitude with exchange 
of a photon is divergent in the limit t — > 0. It should be noted that for large and finite s the 
minimum of \t\ is small but not exactly zero. This fact makes the amplitude finite, it leads to 
that the differential cross-section in the forward direction increases linearly with s and can be 
determined completely. We find that at energies relevant to HERA the contribution from photon 
exchange is actually dominant in the extremely forward direction. This may exclude the possibility 
to observe effect of Odderon exchange in this kinematical region. However, the contribution from 
photon exchange decreases with increasing \t\ more rapidly than that of Odderon exchange, it is 
possible to identify Odderon exchange at \t\ which is not very close to its minimal value. 

Our work is organized as the following: In Sect. 2. we use NRQCD factorization for rj c and 
then formulate the problem of calculating the S-matrix element as the problem of solving wave 
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functions of c quarks in a background field of soft gluons. In Sect. 3. we solve the wave functions by 
expanding them in m~ l . The relevant part of solutions are given. In Sect. 4. we use these solutions 
to derive the ^-matrix element for diffractive photoproduction of J/ty. The results derived in this 
way are exactly the same as those derived by the assumption of two-gluon exchange in a special 
gauge [IB]. In Sect. 5. we derive the results for r\ c production and discuss our results. In Sect. 6. we 
study the contribution through exchange of a photon. Sect. 7. is our summary. 

2. Soft gluon exchange 



We consider the diffractive process: 

l{q) + h(p) -> h(p + A) + Vc (k) (1) 

where the momenta are indicated in the brackets and k = q — A. The Mandelstam variables are 
defined as s = (q + p) 2 and t = A 2 . We consider the kinematic region where \t\ is at order of Aq CD 
and each component of A is at order of Aqcd- The initial photon is real with q 2 = 0, h is any 
light hadron whose mass m is at order of Aqcd- Throughout this work we take nonrelativistic 
normalization for i] c and c-quark. The process undergoes like that the initial photon first splits into 
a cc pair, this pair exchanges soft gluons with the light hadron and forms r\ c after the exchanges. 
A typical diagram for this is given in Fig.l. The S- matrix element can be written as a sum over 
contributions with n-gluon exchange: 



(f\S\i) = (2vr) 4 5 4 (g - A - Jfe) / (Px^xz ■ ■ ■ d A x n d A x J 



d A k\ d A k 2 d A k n d A qi 



(2tt) 4 (2tt) 4 (2tt) 4 (2tt) 4 
.^a 2 ... W2 ...^ ((?) ftj K fe2j . . . k^e-^i^Uxy^m 

■e^ + - ik —(h(p + A)\G^{ Xl )G^(x 2 ) ■ ■■GZ(x n )\h(p)}, (2) 

where i and j stand for color- and Dirac indices for the Dirac fields c{x) and c(0) respectively. 
^aiaz-antniim—nn ^ ^ . . . ^ [ s ^ e amplitude for that the initial photon splits into a cc pair 

and this pair emits n gluons. After the emission the pair is in general not on-shell. In the rest 
frame of r] c the c- or c-quark moves with a small velocity v c . We can expand the matrix element 
(rj c (k)\ci(x)cj{0)\0) in the small v c with NRQCD fields. The leading term in the expansion reads: 



<77 c (/0|Q(x)9(0)iO> 



1 — 7 • v 1 + 7 • v 

n 75 ^ 



J' 



e^ x ( Vc \^ X \0)-(l + O(v 2 )), (3) 



where v = k/M Vc is the velocity of r] c , the matrix element is defined with NRQCD fields in the 
rest frame of T] c , where if}' and x are NRQCD fields of two components, ip^ or x creates a heavy 
quark or a heavy antiquark respectively. The matrix element can be determined from the decay 
width for r] c — > 77: 

r(r/ c -77) = f^l(^ t x|0)| 2 . (4) 

We will neglect higher orders of v c and only keep the leading term in Eq.(3). In this approximation, 
the c- and c quark after emission of n gluons is on-shell because of the projection operators 
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(1 ± 7 • u )/2. They carry the same momentum m c v and the mass of i] c is approximated by 
M Vc = 2m c . Using the spinors u(m c v, si) and v(m c v, s 2 ) one can write the projection operators 
as : 

1 + 7-u or ( w \ 



2 

1 — 7 • v 

= —Ill-U 

2 

S2 



= -m c v°J2v(m c v, si)v(m c v, s 2 ), (5) 



where s 1 and s 2 represents the spin state of c- and c quark respectively. Using Eq.(3) and Eq.(5) 
we obtain: 

(f\S\i) = (2n) 4 5 4 (q - A - k)-^=(m c v°) 2 {r] c \^x\0) £ H™cV, s 2 )>y 5 u(m c v, s ± ) 

d 4 k x d 4 k 2 d 4 k n 



^ 1 f 4 4 4 f df%(f% 

X J -Jdx 1 d x 2 ---dx n j 



(2tt) 4 (2tt) 4 (2tt) 4 
■u(m c v, Sl )A aia2 "' an ' flW2 '" ,ln (q, m c v, k±, k 2 ,--- k n )v(m c v, si) 
•e*"* 1+ "^(/i(p + A)|^(x 1 )(^(x 2 )...G^(x B )|/i(p)>. (6) 

Clearly, the amplitude u(m c v, s 1 )A aia2 "' ar "^ 1 ^ 2 "'^ n (q, m c v, hi, k 2 , ■ ■ ■ k n )v(m c v, si) is just for the pro- 
cess in which the initial photon splits into a cc pair, this pair emits n gluons and becomes on-shell 
when the emission is completed. Hence the ^-matrix element can be written as: 

(f\S\i) = -i{2n) 4 5 4 {q - A - k)—^=(m c v°) 2 (r] c \^x\ty £ v{m c v, s 2 )~i b u(m c v , s x ) 

6Vf u 



Sl,S2 

■(h(p + A), c(m c v, si), c(m c v, s 2 )\J?(0)\h(p))e^ (7) 

where e is the polarization vector of the initial photon, is the charm quark part of the electro 
current, it is J£ = eQ c c^c with Q c = 2/3. With the approximation in Eq.(3) the production 
process can be viewed as a two-step process, in which a on-shell cc pair is produced, then this 
pair is converted into r\ c . The probability amplitude for the conversion is the matrix element 

With the c- and c quark in the final state one can apply the standard LSZ reduction formula 
for the matrix element: 

(h(p + A), c(m c v, si), c(m c v, s 2 )\c(0)^c(0)\h(p)) 

= _L f d 4 xd 4 ye m ^ x+ ^u(m c v, Sl )(i 7 • d x - m c ) 
Z 2 j 

■(h(p+ A)\c(x)c(0)^c(0)c(y)\h(p))(-ij ■ *d y - m c )v(m c v, s 2 ), (8) 

where Z 2 is the renormalization constant of c quark field. We can evaluate the matrix element 
with help of the QCD path-integral, and perform first the integration over c-quark field, while 
gluon fields and other dynamical freedoms will be integrated later. After the integration over 
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c-quark field we obtain: 

(h(p + A), c(m c v, si), c(m c v, s 2 )\c(0)^c(0)\h(p)} 

= f d 4 xd 4 ye m ' v < x+y) u(m c v, Sl )(i~f ■ d x - m c ) 

Z 2 J 

■{h(p + A)\S(x, 0)7^(0, y)|/i(p))(-i 7 • ^ - m>(m c u, s 2 ), (9) 
where S^x, y) is the c quark propagator defined as 

S(x,y) = -(0\c(x)c(y)\0). (10) 

Because the gluon fields are not integrated at the moment and they can be taken as a background. 
Then the propagator describes how the quark propagates under the background of the gluon fields. 
With the propagator we can introduce a wave function for the c- and c quark respectively. We 
define: 

rj} c {x) = u(m c v,si) J d 4 ye imcV ' y (i-fd y -m c )S(y,x), 

Mx) = f d 4 yS{x,y){-i 1 -"d y -m c )v{m c v,s 2 )e im ^\ (11) 
with these wave functions we obtain for the matrix element: 

(h( P + A),c(m c v,s 1 ),c(m c v, S2 )\c(0)^c(0)\h{p)) = ~{h{p + A)|^ c (0) 7 ^c(0)|%)). (12) 

^2 

It should be noted that the quark fields in above manipulations are bar fields, especially, in Eq.(3). 
In Eq.(3) one should use the renormalized fields in both sides to perform the expansion in v c , it will 
give an extra factor of Z 2 . When we substitute Eq.(12) into Eq.(7), this extra factor Z 2 will cancel 
the Z 2 in the denominator of Eq.(12). For gluon fields G M (x) they always appear here as a product 
of g s G^(x), which can be simply replaced with renormalized quantities without introducing any 
extra factor. With Eq.(12) we get the ^-matrix element: 

(f\S\i) = -^^^(g-A-^-^K^^cl^xlO) 

18 Vf 

• v(m c v, s 2 )>y 5 u(m c v, s 1 ){h(p + A)|^ c (0)7 • e^ g (0)|/i(p)). (13) 

Similar results can also be obtained if we consider diffractive photoproduction of J /iff: 

l(q) + h(p) ^h(p + A) + J/iff(k). (14) 
The expansion in the small v c corresponding to Eq.(3) reads: 



(W|Q(x) Cj (y)|0) = J—jm^+y) 

6\/f u 



■(J/^<T-ej(v = 0) X \0) + ---, (15) 
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where e,j is the polarization vector of J/ty. The NQRCD matrix element is related to the decay 
width of J/W -> £+£-: 

Y{J/i, -> e+e") = « e 2 m Qc^KWI^ • e(v = 0) X |0>| 2 . (16) 

With these one can obtain the ^-matrix element for the case with J/ty: 

(f\S\i) = ie(2*)W(q- A-fe)-^=M 2 (J/^|^ ■ e(v = 0) X |0> 

18v v 



• v(m c v, s 2 h ■ e*j(v)u(m c v, s 1 )(h(p + A)|^ c (0)7 ■ e^ e (0)|/i(p)). (17) 



The wave functions defined in Eq.(ll) satisfy the Dirac equations: 

ifj c (x) ■ D + m c ) = 0, 
(ij ■ D — m c ) ipc{x) = (18) 

with the boundary conditions: 

fo(x) - ^m^.^je™ 1 " (19) 

for the time x° — > oo. D M = <9 M + ig s G^{x) is the covariant derivative. These wave functions 
describe the behavior of c- or c quark moving under a background of gluon fields. As discussed 
before, the exchanged gluons between the cc pair and the light hadron h are soft, this results in 
that the background gluon fields are mainly with long-wave lengthes and have a weak dependence 
on x, the dependence is characterized by the scale Aqcd- In the heavy quark limit, m c ^> Aqcd, 
this enables us to solve the equations by expanding the wave functions in m~ l . Using the solutions 
we can obtain the matrix element in Eq.(12) in terms of gluon fields. Then we integrate over the 
gluon fields, i.e., complete the QCD path integral mentioned after Eq.(8), and finally we can 
get the S'-matrix element for the diffractive process. It should be noted that the propagator in 
Eq.(10) is a Feynman propagator. The wave-functions defined in Eq.(ll) do satisfy the Dirac 
equation with a background field of gluons, but they do not have a simple boundary condition like 
in Eq.(19). However, it was shown [T3] that the Feynman propagator in Eq.(10) can be replaced 
with a advanced- or retarded propagator, if the background field varies enough slowly with the 
space-time. Therefore one can have solutions of Eq.(18) and they satisfy the boundary conditions 
in Eq.(19). In our case we will make the expansion in m~ l and in this expansion the c-quark and 
c quark are decoupled, and it results in that the Feynman propagator will automatically be an 
advanced propagator for c and c in our case. 

3. The wave functions in the expansion in m^ 1 

In this section we will solve the Dirac equations up to order of m~ 2 for our purpose. With the 
invention of HQET the expansion in m~ l is now quite standard. Taking tp 5 as an example, we 
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can decompose the wave function as: 

fc(x) = e im ^ x (h + (x) + h-(x)) , (20) 

where 

7 • vh + (x) — h+(x), 7 • vh-(x) — — h-(x). (21) 

To present the solution for ip 5 we introduce some notations. For any vector A one can decompose 
it as: 

A» = vAv ll + A£, with v ■ A T = 0. (22) 
Using the Dirac equation for ip 5 one can express h+(x) in term of h-(x) as a series of m~ 



-n. 

"c ' 



h+(x) = - -i-f-D T h-(x) 

2m c — iv ■ D 



1 ^ r -D T h (,•) • 0{\) (23) 



2m c Ami J m c 



while h-(x) is constrained by the equation: 



c / ) c 



This equation can be solved by expanding h-(x) in to" 1 , and eventually we obtain the wave 
function. It is straightforward to solve this equation. For convenience we define a gauge link V(x) 
with = ujv^ + x% as: 



V(x) = Pexp 



roo 

-ig s / drv ■ G(r,x T ) 

J U) 



(25) 



where we denote the x-dependence of G M (x) as G^(uj,x t ) and P stands for path ordering. The 
leading order results for wave functions read: 

ip e (x) = e im ^ x V^(x)v(m c v,s 2 ) + 0(—), 

ui c 

$J X ) = e imcV - x u(m c v, Sl )V(x) + 0( — ). (26) 

m c 

With these wave functions at the order of m° c it is easy to find out that the ^-matrix element in 
Eq.(13) or Eq.(17) is zero because VV^ = 1. To go beyond the leading order we make a gauge 
transformation: 

ip e (x) -> V(x)ij) B (x), 
$ c (x) -> V f (x)^ c (x), 

GHx) -> Wx^fx^fx) - — Vfx^ytfa;). (27) 

The S 1 - matrix element is invariant under the gauge transformation. The transformed gauge field 
has no component along the direction v, i.e., v ■ G(x) = 0. This is equivalent by taking the gauge 
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v ■ G(x) = 0. To avoid introduction of too many notations we will use the same notations for 
transformed fields. The fields below should be understood as transformed fields or fields with the 
gauge v ■ G(x) = 0. Now for gluon fields we have: 

v ■ G(x) = 0, VpG^ix) = v^G v {x). (28) 

At the order of m~ l the solutions for the wave functions read: 



c 



1 + _ 7 • D T - — I dr {l ■ D T f) V (m cV , s 2 ) + O(^) 

= e^/i + ^L r d rG(r,x T )-G(r,x T ) + --]v(m c v,s 2 ) + 0(^) 
{ 2m c Juj J m% 

= e m ^u(m c v,s 1 ){l--^ 1 -<D T --^- [°° dr^ ■ JD T ) 2 } + O(^) 
I 2m c 2m c Juj J m l c 

= e m ^ x u{m c v,s 1 )\l + ^- [°° drG(r,x T ) ■ G(r,x T ) + ■■■} + 0{\), (29) 
[ 2m c Juj J m z c 

where in the {• • -}'s with • • • we have detailed terms, which will lead to contributions to the S- 
matrix element in the case with J/^f, and • • • denotes irrelevant terms. At this order the S-matrix 
element in the case with r] c is zero. To obtain the nonzero S'-matrix element one needs the solutions 
at order of m~ 2 . 

At order of m~ 2 there are many terms for the wave function ipc{%)- However, only one term 
will lead to a nonzero contribution to the S'-matrix element with r) c . We will only give this term 
in detail for the solutions, other irrelevant terms are denoted by • • •. For ipc{%) we have: 

tcix) = jl + -L {•••}- -L QH dri ( 7 . D T f jf° rfr 2 ( 7 • D T f + ■ ■ •) j v(m c v, s 2 ) 

+0(— 3 ), (30) 

one can use the identity 

( 7 • D f = D ■ D + (31) 

to indicate the relevant term more clearly. The part of wave functions relevant to the case with 
r] c reads: 

i>- c {x) = e m ^ x {l + — {...}- JL[[ °° dn 1°° dr 2 {g s G-G-\a, v G^){n,x T ) 

77l c 4:Tfl c Juj J t\ £ 

■{g s G ■ G - \a, v Gn{r 2 , x T ) + ■ ■ -}}v(m c v, s 2 ) + 0{\), 
$ c {x) = e ^-u{m c v,s l ){\ + — {•••}- f°° dn l°° dr 2 {g s G ■ G - l -a^)(r 2 ,x T ) 

Ul r 4?71_ Juj J ri 2 



■{g s G ■ G - i^G^Xri, xt) + •••]} + O(^). (32) 



c 



With the terms for the solutions, given in detail in Eq.(29) and Eq.(32), we can evaluate the 
S- matrix element for J/ty and r\ c . 

4. The S- matrix element for J/ty 

In this section we use the solutions of wave functions to derive the ^-matrix element for 
diffractive photoproduction of J/^f. Using Eq.(29) we obtain the term in Eq.(17) 

(m c v ) 2 J2 v(m c v, s 2 )7 • e*j(v)u(m c v, Si)(/i(p + A) ^(0)7 • eip 5 (0)\h(p)} 

Sl,S 2 

= l J^e* J .e f°° du(h(p + A)\G a >»(u,0)G;(u,0)\h(p)). (33) 
Using translation covariance and 

/ due™* = (34) 

Jo v-A + iO+ K J 

where + is a positive infinitesimal number, we obtain 

y ) 2 v(rn c v, s 2 ) 7 • e*j{v)u{m c v, si)(% + A)|^ c (0) 7 • etpc(0)\h(p)) 



(m r v ^ 2 



S1,S2 

,2 



= ~— 2 e* ■ e(h(p + A)|G^(0, 0)G£(0, 0)|%)>- (35) 
m c 

In the above equation we has used the kinematics for the diffractive process: 

v ■ A = -m c (l + » -m c . (36) 

For the gauge field with v ■ G = as in Eq.(28) one can relate it to the gluon field strength: 

/oo 
dr9(T-cu)v u G^(T,x T ), (37) 
-00 

where the step function is defined as: 

too dw fe-*" _ [ l,t>0 

m ~ LootoZTfW ~ \ o,t<o (38) 



With this relation one can show that 

, (0.0)\h(n)\ = — f dz- 

+ z~iO+ l-z-iO+ 



(Mp + a)|G7(o,o )G ^(o,o)|Mp)) = ^J^ TT h W ■ T-hi ( 

/rjf 
—gy mcZT vMKv + A)\G a ^(r, 0)G a ' p (-r, 0)\h(p)). (39) 



Finally we obtain the S-matrix element for J/\&: 

(f\S\i) = ^(2vr) 4 5 4 (g-A-fc)-^(J/^|^ CT . £ ( v 

r 1 1 

• / <i,2 ■ ■ Fr(z), 

J l + z-iO+ l-z-iO+ K J 

with 

/d,T 
— g y m ^ v v u ( h (p + A)\G a ^(Tv)G a '%(-Tv)\h(p)). (41) 

This result is exactly the same as that given in [IBj, where we use perturbative QCD by taking 
only two gluon-exchange in the gauge v ■ G = 0. Here we derive the same result without using 
perturbative QCD, instead we only use the expansion in m~ l in an arbitrary gauge. The above 
results are expressed with the transformed gauge field, the gauge transformation can be found in 
Eq.(27). If we express the results with the untransformed gauge field, a gauge link will appear 
automatically between the two operators of gluon strength fields. The gauge link is along the 
direction of v with the gauge fields in the adjoint representation, it starts at x — —tv and ends at 
x = tv. With the gauge link the results are gauge- invariant. It is interesting to note that models 
with two-gluon exchange are widely used for diffractive photoproduction of vector meson. Our 
result here shows that such a model is a correct approximation in the heavy quark limit. 

In our approach, although a formal expansion in m~ l is employed, but the true expansion 
parameter is (m c v ■ k) , where k is the momentum of a exchanged gluon, this can be realized 
by inspecting the lagrangian of HQET and it is discussed in detail in [IB] . In this expansion, 
transversal momenta of exchanged gluons are neglected at the leading order. Hence the exchanged 
gluons will not resolve the structure of a heavy quarkonium in the directions transversal to v. The 
situation here is similar to the multi-pole expansion for gluon fields in hadronic transitions like 
V>' -> J/"® + 7T + 7T in [T3] . 

For sufficiently large beam energies, the dominant contribution to the correlation function 
Fr{z) is from the standard gluon-operators with twist 2 and these operators are used to define the 
gluon distribution of h. If we only take the dominant contribution, then the 5*-matrix element is 
related to the generalized gluon distribution. However, one can show that the S'-matrix element 
can be related to the usual gluon distribution gh(x) at x — 2m 2 c / s for t — > 0, this is different than 
that in previous approaches [20 , details can be found in [IB] . 

In [18J it shows that the predicted cross-section has large deviation from experimentally mea- 
sured for J/ 1 ®, while for T the agreement is fairly good. One of possible reasons for the large 
deviation can be corrections from higher orders in m" 1 , because m c may not be large enough, while 
for T the 6-quark mass is heavy enough to have reliable results from leading order. At leading 
order the exchange is of two gluons as indicated in Eq.(29). It is interesting to note that with the 
method developed here, it is possible to resume many exchanges of this type of two gluons, the 
resummation will reduce the corrections from higher order in m~ l . Works along this direction are 
under progress. 

5. The ^-matrix element for rj c 



o)x|o) 



(40) 
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Using the wave functions given in Eq.(30), it is straightforward to obtain the S'-matrix element 
for i] c . We obtain the corresponding part given in Eq.(33): 

(m c v ) 2 v(m c v, s 2 )~f 5 u(m c v, Si){h(p + A) ^(0)7 • eij) B (0)\h(p)) 

S1,S2 

,3 



9 s 



/oo /*oo 
dn / dT 2 6(T 1 )6(T 2 - n) 
-00 J — 00 



■d abc (h( P + A)\ {G a ^{nv)G h p {nv)G c ^ \t 2 v) + (n - t 2 ,t 2 - n)} |%)>. (42) 

With help of Eq.(37) one can transform the gluon fields G M into the field strength G^ . After 
some algebra manipulation we can express the result as: 

(f\S\i) = ie(2nf5\q -A- k) 9s (^^K^e'v? [ du x du 2 ■ T^{^ 2 ) 

(2m 6 c yv [i J 

{cj! + iO+ cu 2 + i0+} { m c + u x + uj 2 - i0+ uj x + uj 2 + i0+ } ' ^ 

with 

^(uj 1 ,lo 2 )= f pL^le-i^-iu,^ . d abc^ p + A)|G' a ' CTQ (r 1 ^)G'^ Q (r 2 ^)G' c '^(0)|/i(p))^ p . (44) 

Again, in the definition of T^ioj^ oo 2 ) the gluon fields are transformed fields in Eq.(27) or those in 
the gauge v ■ G(x) = 0. If we use the untransformed to express Eq. (44), the gauge links discussed 
in the last section will automatically appear between the field strength operators. Therefore, the 
^-matrix is gauge invariant. From this result one can also see that in the heavy quark limit, 
the exchange of three gluons in the gauge v ■ G(x) = is responsible for the process and the 
polarizations of the three gluons are all transversal to the direction v, and these three gluons 
are emitted by c- or c quark. In other gauge, because the appearance of the above mentioned 
gauge links, there are additionally exchanges of gluons with infinite numbers, those gluons have 
polarizations proportional to the direction of v and their effects are included in the gauge link. 

If the beam energies are large, which is really required by that t ~ ^qcdi the vector v 
approaches to a light-cone vector, the dominant contribution for T^ v can be identified. For 
convenience we take a coordinate system in which the photon moves in the — z-direction and the 
hadron h in the z-direction. We introduce a light-cone coordinate system, components of a vector 
A in this coordinate system are related to those in the usual coordinate system as 

^ = {A\ A-, A T ) = A\ A^ . (45) 

We introduce two light-cone vectors n = (0, 1, 0, 0) and I = (1, 0, 0, 0) with I -n — 1. The momenta 
in the process can be approximated in the limit s ~ 2k~p + — > 00 as: 

k" = (0,AT,0 T ) = fc"n", 

9 

777 

p" = (p + ,— ,0 T )«p + P, 
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(k-AY fm c k~ n \ k~ u 
v = \, ~ M,^— ,0t ^— n\ (46) 



Mj/^ \k~ 2m c J 2m c 

with the above approximated momenta and re-arrangement of variables the dominant contribution 
of the 5-matrix element reads: 

(f\S\i) « ie(2n) 4 5 4 (q- A - k) ^-^(rj c ^ x \0)e ap ^n p F [ d^du 2 ■ W> 1; u 2 ) 

36m^Vf J 

(47) 



. u\ + i0 + uo 2 + i0+ J 1 2x c + uoi + uo 2 — i0 + tui + co 2 + «0 + 
with x c = 2ml/ s an d 



W>i,w 2 ) = — / ■7r-7r-e" la ' iriP W2T2P • d (h(p + A)\G a,+a (T 1 n)G ,+ a (r 2 n)G c,+,/ (0)\h(p)) , 

(48) 

where a>i and co>2 is dimensionless and i/ can only be 1 or 2. The function W(u>i,uj 2 ) with u — 1,2 
is invariant under a Lorentz boost along the z-direction. W is defined by the twist-3 operator in 
the matrix element and its dimension is 1 in mass, hence it is proportional to A, where A can be 
one of the small scales like Aqqd, \f—i, etc.. This is our main result for any light hadron. If one 
takes the approach with exchange of a Odderon for the process, the effect of the Odderon is then 
represented by W\ which is defined in the framework of QCD. 

Now we consider the case that h is a proton, for this case W ! (i = 1,2) can be parameterized 
with four functions in general: 



W = \u(p + A) 
p+ 



f ■ +i, f p A * , f A V . , p + Y 

+h — 2~ + / 3 + 

m m m 



u{p), (49) 



where u and u is the spinor of the proton, m is the proton mass. A similar decomposition of quark 
operator can be found in [2U E2j, a general discussion about how to write these form factors in 
similar cases is given in [21]. The functions fi(i = 1, 2, 3, 4) are with dimension 1 in mass and are 
with variables: u>i, uj 2 , t and £ = — A + /(2p + + A + ). These functions are proportional to A. The 
differential cross section can be expressed with these four functions, the expression is too long to 
present here. However one can look at how the differential cross-section depends on dimensional 
quantities, like s, t, m c etc. . We have for t — ► 0: 



-rriVc) ~ a, 3 2 • —f{s, t, m c ), (50) 

dt m± mi s z 



the dimensionless function / is determined by the four functions given above. It is interesting to 
note that in the limit A 1 ' — > da/dt does not approach to zero. In this limit the contributions 
from f 2 and indeed vanish , but the contributions from fx and f± survive in the limit. This 
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fact is in contrast to the results from EH E] ■ This result can be understand as the following: 
Because the initial photon has the helicity A 7 = ±1 and r] c has the helicity 0, the proton will 
change its J z with ±1, the component in ^-direction of the total angular momentum. This change 
can be made with change of l z , the component of the orbital angular momentum, whose effects 
are parameterized by f'2 and f'3 and the change is represented by A 1 in Eq.(49). It is clear that 
in the limit A — > the contributions from f'2 and ^3 go to zero. However, this change can be 
made by changing the helicity of the proton, which can be identified in massless limit as s z , the 
component in z-direction of the proton spin, effect due to this is given by f\, the term with f± 
conserves helicity, but s z can not be exactly identified with the finite mass m as helicity, if m is 
exactly zero, the term does not exist in the limit A M = 0, hence the nonzero contribution from 
fi in the limit A M = is due to the finite mass m. These contributions survive in the limit and 
there are in general no reasons why they should vanish. If one uses perturbative QCD and takes 
a set of light-cone wave functions of the proton to calculate W or the impact factor in EI] for 
interactions between three gluons and the proton, where one takes a proton as a bound state of 
the three quarks uud and the three quarks are in s-wave, then one would find that f\ = because 
perturbative QCD conserves helicities. However, if one realizes that fact that the three quarks can 
have orbital motions, whose effects are parameterized with wave- functions given in [22], then fi 
is not zero. In general perturbative QCD is questionable for this calculation, because the gluons 
are soft. However, this indicates f\ 7^ 0. 



If one thinks that the gluons, whose effect is represented by W, may be generated by quarks 
inside the proton, one may relate W l to various quark distributions, in which one avoids to use 
wave functions of the proton. Because W % is defined with a twist-3 operator and contributions 
from quark operators of twist-4 are presumably suppressed, then these gluons are only generated 
by one quark. This can be illustrated with Fig.2, where we use perturbative QCD as a guide and 
only three gluon lines are drawn. In Fig.2. the black dot represents the gluonic operator. Up 
to twist-3, all gluon lines can only be attached to a single quark line, if gluon lines are attached 
different quark lines, it will results in contributions at twist-4 level or higher. Then W 4 will be 



+■ 




h{p) 



h{p + A) 



Fig.2 



Figure 2: Gluons are generated by a single quark in a proton, see text 
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proportional to the quark density matrix: 



(pip + A)|gi(-rn)^(— irn)|p(p)), (51) 

where i, j stand for Dirac indices. If the calculation can be done in this way with perturbative 
QCD, then one will find that W ! are related to various nondiagonal quark distributions, which 
may be found in [2UI22JI2]- ^ is interesting to note that the helicity flip quark distributions will 
contribute to f\ in this calculation. 

It is constructive to look at how the differential cross-section of J/^f depends on dimensional 
quantities to compare with that of rj c in Eq.(50). From JH] we have for t — > 0: 

— (W)~a fl -j 2 -\9h{x c )\ , (52) 

in this case we know that the dimensionless function corresponding to / in Eq.(50) is determined 
by the usual gluon distribution with x c defined after Eq.(47). For x c — > the gluon distribution 
behaves like x~( 1+ ^ with /? > 0, this results in that the cross-section for J/ty increases with 
increasing s. If we assume that / has a similar dependence on s as gh{x c ), we can roughly obtain 
the estimation: 

d&(„ \ \/ m L/,t,,ln\l2 a2 



dt 



(Vc) \(Vc\^x\0)\ 2 A 



- a s — — ; — : — t : : — ; — ttt: • — ;r = 0.008 (53) 

§(J/V>) s |(J/^|^t o -.s(v = 0)x|0)| 2 ml 1 J 

where the NRQCD matrix elements are determined by the decay rj c — ► 77 and J/^> — >• e + e~, 
respectively, the small scale A is taken as 300MeV, other parameters are taken as a s = a s (m c ) ~ 
0.3 and m c = 1.5GeV. The estimation is rough but it shows that the production r\ c is suppressed 
in comparison with J/ty. Using the experimental data from HERA for J/?/>[26 we estimate for 

t -> or y/\t\ ~ A = 300MeV: 

^(r /c )~2nb/GeV 2 (54) 
dt 

at y/s = lOOGeV. This rough estimation gives a larger number than those in jUJ CHI ITTj . For 
example, in jTU] da/dt is about O.lnb/GeV 2 for yj\t\ ~ 300MeV. In jlj it is pointed out that the 
proton-Odderon impact factor used in jHl E5] may be overestimated, this will result in that the 
prediction for da/dt in j^JCni may also be overestimated. This makes that the difference between 
our estimation and that in CHI becomes larger than mentioned above. It should be noted that 
our results in Eq.(43) and Eq.(47) can have large uncertainties. Because c-quark may not be 
heavy enough, corrections, like relativistic correction and those from higher orders in m" 1 , can be 
significant, as already discussed the case for J/ty. For bottomonia, those corrections are expected 
to be small because of the large 6-quark mass. 

In our work we show that at the leading order of m~ l the dominant contribution is from four 
distribution functions defined with gluonic operator at twist 3. At higher orders operators at higher 
twist will appear, because there are some higher-twist effect in the 3-gluon change and more than 3 
gluons can be exchanged. One can use some resummation techniques for gluon exchanges instead 
of the expansion m" 1 , like by using BFKL equation. In our case with Odderon the equation for the 
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resummation is the Bartels-Kwiecinski-Praszalowicz (BKP) equation When the resummation 
is used, the effect of higher twist contributions in the 3-gluon exchange, including that from twist 
3 contribution, will be taken into account. At moment, no information is available for the matrix 
element with the twist-3 operator in Eq.(48), hence a detailed prediction is not possible. 



6. Contribution from Photon Exchange 



In this section we study the contribution from photon exchange instead of exchange of many 
gluons in Fig.l. The contribution is represented by Fig.3. We will only concentrate on the 
differential cross-section in the extremely forward direction, i.e., in the limit t — > 0. In this 
limit, the contribution is divergent at first look because the exchanged photon becomes soft. 
But, it should be noted that for large and finite s, t can never be zero exactly, for large s, i.e., 
s ^> M 2 c « (2m c ) 2 , t approaches to —\t\ m i n in the nearly forward direction: 

m 2 (2m r ) 2 , . 

f ~ —\t\ \t\ ■ Ri f55^ 

L \ L \mim \ b \min ~ 5 V""/ 

S 

this minimum value will acts as an infrared cutoff and makes the S"-matrix element finite in the 
forward direction. Exchange of soft gluons can happen in combination with the photon exchange 
in Fig.3, it will lead to contributions in the forward region which are suppressed by A 2 /m 2 and 
these contributions can be neglected. 




Vc(k) 



hip + A) 
Fig.3 

Figure 3: One of the two diagrams for photon exchange 
The calculation of the S-matrix element is straightforward, we obtain: 

(f\S\i) em = z2(eQ c ) 2 (27r)V(g- A-A;)^(^I^X|0)W^A^ 

-^■(h{p + A)\no)\h(p)), (56) 



A 2 + 2m c v ■ A A 

where J M is the operator of electric current, Q c is the charge fraction of c-quark in unit e. For 
t — ► —\t\ m i n , one may use the method developed here to define a wave function for the c- or c 
quark propagating in a background field of soft photon to calculate the S- matrix element, the same 
result can be obtained if one neglects A 2 in comparison to 2m c v ■ A in the above equation. In this 
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way one can also show that the contribution from a photon exchange combined with exchanging 
soft gluons is suppressed by A 2 /m 2 c . For h being a proton, the matrix element is decomposed with 
Dirac- and Pauli form factor: 

{hip + A)\J»(0)\h(p)) = uip + A) | 7 ^i(A 2 ) + io^^SJ u(p), (57) 

for t = —\t\ m i n the differential cross-section can be obtained directly. The dominant contribution 
reads: 



da „ <^Qt J_ F 2 

dt lt ~ ltlmin ~ 8m! m? m 2 



where we use F\ m e for the small momentum transfer. For large t, because the form factor F\ (t) 
falls like t _2 |2ZI and the form factor ^(t) falls like t~ 3 |28j. the differential cross-section will fall like 
t~ 8 . A model calculation in |U] for gluon exchange predicts that the differential cross-section falls 
t~ 2 . Hence, at large t the differential cross-section with the photon exchange is highly suppressed. 
Because the detailed information of these form factors are not available, numerical result of this 
contribution to the total cross-section can not predicted. For small t, the differential cross-section 
increases as at t — — \t\ m i n it increases with s linearly, this is because the exchange photon 
becomes soft and its invariant mass decreases with s~ l . From Eq.(58) the differential cross-section 
in the forward region can be large, although it is suppressed by o? em - The numerical value of da/dt 
can be worked out, we have: 



da , / \fs \ nb . . 

8 4X10 (l&v) CSV*" <59) 

where we used m c = 1.5GeV. With this estimation one can see that at = lOOGeV, the 
contribution from photon exchange is already at the same order as the estimation given in Eq.(54). 
It is also larger than the numerical values given in jUJ E3 EJ. Hence, if the production of r] c is 
observed in the extremely forward region the contribution from photon exchange is dominant for 
the production at energies relevant to HERA. This brings the difficulty to identify the effect due 
to Odderon exchange in the extremely forward region. For s larger than those at HERA the 
contribution will still be dominant, if the contribution from Odderon exchange increases with s 
as s a for a < 1. Since the contribution from photon exchange decreases with increasing \t\ more 
rapidly than that of Odderon exchange, it is possible to identify the effect of Odderon exchange 
for suitable \t\ which is not very close to |t| m in- It requires a detailed study to identify kinematical 
regions for hunting Odderon. 

7. Summary 



In this work we have studied diffractive photoproduction of rj c . Taking charm quark as a heavy 
quark, the nonpertubative effect related to r] c is represented by a NRQCD matrix element and 
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r] c can be taken as a bound state of cc quark. Then the production can be imagined as that 
the initial photon splits into a cc pair, this pair forms the r] c after exchange of soft gluons with 
the initial hadron in the forward region. The problem of exchange of soft gluons can be studied 
with HQET and a systematic expansion in m~ l can be employed. We find that the ^-matrix 
element can be expressed in terms of wave-functions of c- and c quark, which propagates under 
a background field of gluons. The background field is dominated by components with long-wave 
lengthes corresponding to that the exchanged gluons are soft. By solving these wave functions 
with the expansion in m" 1 , we can express the result for the S- matrix element with quantities, 
which are well defined in QCD. The S'-matrix element in the forward region is expressed in the 
case of proton with four functions, which are defined with a twist-3 operator of gluons. If one 
thinks that the exchange of Odderon is responsible for the production, then the effect of Odderon 
is represented by these functions. For completeness, we also derive the S'-matrix element for 
diffr active photoproduction of J/^, the result is exactly the same as that derived in JH] with an 
assumption. 

Since these four functions are unknown, numerical predictions can not be made in this work. 
However, an estimation in order of magnitude can be given in comparison with the production 
of J/ty . The estimation indicates that the differential cross-section of r] c in the forward region 
is really small, in qualitative agreement with previous studiesjHl HDl IH], where some models are 
used to make numerical predictions. However, from our results the differential cross-section of r] c 
is not zero in the limit t — ^ 0, this is different than those predicted in EH 111 j - The reason for 
this nonzero is that the initial hadron can change its helicity in this limit. 

Instead of exchanging many soft gluons, a photon exchange can also contributes to the produc- 
tion of r] c . In the forward direction, i.e., t — > 0, one may find that this contribution is divergent. 
But, we should note that t can approach to for s — * oo, for large and finite s, the minimum value 
of \t\ is never zero because all particle in the final state are with finite masses. This minimum 
value makes the contribution finite. At this minimum value of \t\ the differential cross-section due 
to photon exchange can be predicted. We find that this contribution is large and is dominant 
for production of r] c . Hence, if 7] c is produced in the extremely forward region in experiment, one 
can not conclude that the r\ c is produced through Odderon exchange and the effect of Odderon is 
observed. At large \t\, it can be shown that the differential cross-section with the photon exchange 
falls like t -8 and can be neglected. Therefore, if rj c is produced with \t\ which is not very close to 
its minimal value, one can identify the production as an effect of Odderon to confirm the existence 
of Odderon in this process. 

Our results for scattering amplitudes are expressed with quantities like NRQCD matrix element 
and matrix element of the twist-3 operator of gluons, they are well defined in the framework of 
QCD and are universal, i.e., they do not depend on specific processes. Nonperturbative methods 
like lattice QCD or sum rule method may be used to determine them. The matrix element of the 
twist-3 operator is unfortunately unknown at the moment, this fact prevents us to make numerical 
predictions. But, with our results one can build more realistic models to have a reliable prediction 
for diffractive photoproduction of r\ c . 
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